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ABSTRACT 

Doppler surveys have shown that the occurrence rate of Jupiter-mass planets appears to increase 
as a function of stellar mass. However, this result depends on the ability to accurately measure the 
masses of evolved stars. Recently, Lloyd (2011) called into question the masses of subgiant stars 
targeted by Doppler surveys. He argues that very few observable subgiants have masses greater than 
1.5M Q , and that most of them have masses in the range 1.0-1.2 M©. To investigate this claim, we 
use Galactic stellar population models to generate an all-sky distribution of stars. We incorporate the 
effects that make massive subgiants less numerous, such as the initial mass function and differences 
in stellar evolution timescales. We find that these effects lead to neglibily small systematic errors in 
stellar mass estimates, in contrast to the 50% errors predicted by Lloyd. Additionally, our simulated 
target sample does in fact include a significant fraction of stars in excess of 1.5 Mq. The inclusion of 
an apparent magnitude limit results in a Malmquist-like bias toward more massive stars, in contrast to 
the volume-limited simulated sample of Lloyd. The magnitude limit shifts the mean of our simulated 
distribution toward higher masses and results in a relatively smaller number of evolved stars with 
masses in the range 1.0-1.2 Mq. Thus, many of the subgiants targeted by Doppler surveys are likely 
the evolved counterparts of massive main-sequence stars. 
Subject headings: stars: fundamental parameters — stars:subgiants 



1. INTRODUCTION 

Studies of the relationships between exoplanets and 
their host stars provide valuable clues about how plan- 
ets form, and also point the way to new discoveries. 
For example, the well-established relationship between 
the occurre nce rate of gas giant planets and host-star 
metallicity dSantos et a ll 120041: IFi scher fc ValentH 120051: 
Uohnson & Appsll2 009) may be an indication that the for- 
mation timescale for close-in giant planets (a < 5 AU) is 
shortened by the metal-enhancement, and h ence dust- 
enha ncement, of protoplanetary disks (e.g. Ilda fc Lml 
2004). For this reason, certain Doppler surveys have 
biased their target lists toward metal-rich stars, which 
has resulted in t he discovery of many of the known hot 
Jupiter systems ([Fischer et al.ll2005t (Bouchv et aLll2005t 
ISato et al.ll2005h . 

More recent Doppler surveys have discovered that 
stellar mass is another key pred i ctor of giant planet 
occurrence (| Johnson et all l2007al l2010al) . This rela- 
tionship is based on Doppler surveys of M dwarfs on 
one si de of the stellar mass range (e.g. Uohnson et al.l 
I2010cl ). and the evolved counterparts of F- and A- 
type stars on the more massive end (iJohnson et al.l 
120071* ILovis fc Mavon [20071: ISato et al.1 12007D . These 
so-called "retired A-stars" exhibit dramatically slower 
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rotation velo cities (Vrotsini) than their main-sequence 
proge nitors ()Grav fc Nagarl fi"985l: Ido Nascimento et al.l 
2000), making them better targets for Doppler-based 
planet surveys compar ed to their F- and A-type main- 
sequence counterparts (|Hatzes et al.l 120031 : IFischer et al.l 
[200l iGalland et al.ll2005r ). 

However, the mass estimates of subgiants targeted by 
Do ppler survey s have recently been called into question 
by iLlovdl ((201 lL hereafter Lll). In an attempt to study 
the effects of star-planet tidal interactions in planetary 
systems with evolved host stars, Lll investigated the ex- 
pected mass distribution of evolved stars near the sub- 
giant branch. By using stellar evolution model grids, 
assumptions about the metallicity distribution in the 
Galaxy, and the form of the stellar initial mass func- 
tion (IMF), Lll concluded that most bright subgiants 
are not the evolved brethren of A-type stars, but rather 
the evolved counterparts of Sun-like stars. This is be- 
cause massive stars evolve much more quickly along the 
subgiant branch than do less massive stars. Lll predict 
that this differential evolution rate, together with the 
IMF, should result in a very small number of subgiants 
with M > 1.5 M in Doppler survey^. 

In this Letter we reassess the conclusions of Lll using a 
simple Bayesian frame work with inp u t from the Galactic 
population models of iGirardi et al.l (j2005j ). As we will 
show, neglecting to account for the different evolution 
rates of stars with different masses results in a negligibly 
small bias in the mass measurements for individual sub- 
giant stars. We als o demonstrate that th e target selec- 
tion criteria of the ([Johnson et al.ll2"010b| ) Keck Doppler 
survey can be expected to result in a higher fraction of 
massive stars than suggested by Lll. 

5 While Lll discuss stellar rotation in great detail, it is this 
evolution rate feature that is his key argument that subgiant masses 
are incorrect. 
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2. ESTIMATING THE MASSES OF SINGLE STARS 

The heart of the problem is that measuring the masses 
of single stars is necessarily a model-dependent proce- 
dure. The most common method of estimating masses 
is to interpolate theoretical stellar evolution grids at the 
positions of various measured stellar properties. Typ- 
ically, the set of parameters used are the stellar effec- 
tive temperature (T e f[), luminosity (L; or absolute mag- 
nitude A4y) and metallicity ([Fe/H]) based on LTE at- 
mospheric inodjls_Jitted_to high-r e solution stellar spec - 
tra (e.g. iValenti fc Fischen [20051: iTakeda et all [2008h . 
The top panel of Figure Q] illustrates a simplified in- 
terpolation in which measurements of T e g and L at a 
fixed [Fe/H] = (red circle and error bars) are com- 
pared to mass tracks from the Yale Rotational Evolu- 
tion Code models (YR EC; solid lines; ITakeda et alj|2507t 
iDemarque et al.ll2008l solid black lines). We also demon- 
strate the effect of metallicity by showing two solar- 
mass tracks with [Fe/H] = {—0.16,-0.50}, which rep- 
resent roughly 2.3-c and 8.0-er deviations from the mea- 
sured metallicity ([Fc/H] = —0.02), respectively. This 
illustrates the large systematic measurement errors that 
would be necessary to mistake a 1.5 M subgiant for a 
less massive star. 

2.1. A Probabilistic Framework 

The probability of a star's mass, M, given its spectro- 
scopic parameters and the selection criteria of a survey 
is given by Bayes' theorem 



P(M\T eS , [Fe/B], L, M v , B -V,l) cc 
P(T eff , [Fe/H], M v | M)x 

P(M | M V ,B-V,T) (1) 

The left-hand side of the proportionality is an expression 
for the posterior probability distribution of the stellar 
mass given a spectroscopic estimate of the stellar effec- 
tive temperature T c g, metallicity [Fe/H], and bolometric 
luminosity, L. In addition to the spectroscopic parame- 
ters, we also have additional information I, which in our 
an alysis is provided by the galactic population models 
of iGirardi et al.l ([2005). The term I encodes informa- 
tion about the stellar IMF, stellar evolution models, and 
the distribution of ages and metallicities a s a fu nction 
of Galactic scale height (see iDawson et al.l I2012L for a 
similar application). 

The right-hand side of Eqn. [T] is the product of two 
probabilities. The first is the likelihood, which relates 
the probability of measuring the spectroscopic properties 
of the star given various choices of the stellar mass from 
stellar evolution models. The second term describes our 
prior knowledge about the distribution of stellar masses 
for stars throughout the Galaxy with a given range of 
photometric properties. 

It is common for investigators using model grid in- 
terpolations to focus solely on the likelihood term, be- 
cause the maximization of the likelihood is directly re- 
lated to the concept of "chi-squared minimization" when 
the measured parameters are normally distributed. This 
can be seen by taking the logarithm of the likelihood, C, 
with normally-distributed measurement uncertainties on 
the spectroscopic parameters in Eqn. [T] 
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Fig. 1. — An illustration of the interpolation of model-grids at 
the position of a subgiant's spcctroscopically-measured T c g and L. 
Top: The subgiant branch of the theoretical Hertzsprung-Russell 
diagram near the base of the red giant branch. The solid lines 
show the evolutionary paths of stars of various masses. The red dot 
shows the position of a specific subgiant and the error bars show the 
68.2% confidence region of its spectroscopic parameters. The dot- 
dashed lines show two other metallicities for the Solar-mass track, 
corresponding to 3 and 10 times larger than the typical errors in 
spectroscopically measured metallicities. The dashed box traces 
the region shown in the lower panel. Bottom: A zoomed-in region 
around the subgiant's measurements, with the mass tracks shown 
sampled at a uniform time-spacing of At = 1 Myr illustrating the 
different evolution rates for subgiants of various masses. 



C = In [P(T eff , [Fe/H], L)} 

= C- ^(XT cff + Xfpe/H] + Xl) 



where, e.g. 



T c ff(M)-r off ,, 

CTeff 



(2) 



(3) 



Minimization of the \ 2 terms maximizes the likelihood 
of the measurements as a function of M. However, this 
least-squares approach neglects prior information about 
the distribution of stellar masses. 

2.2. The Stellar Mass Prior 

Even though the measurement errors for stellar prop- 
erties may be symmetrically distributed across several 
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Fig. 2. — A typical mass probability distribution for a 1.7 Mq 
subgiant based on atmospheric parameters interpolated onto stellar 
evolution grids is (likelihood; solid red line), and the prior distribu- 
tion for stars with comparable colors and magnitudes (solid black 
line). The resulting mass posterior is the product of the likelihood 
and prior (dashed blue line) and has a mean of 1.65 Mq (blue ar- 
row) , which is roughly 3% lower than the initial estimate without a 
prior (red arrow). The scaling of the distributions does not matter 
since such factors scale out of the product and do not affect the 
centroid or symmetry of the final posterior distribution, hence the 
proportionality in Eqn. [T] 

mass tracks, there is not an equal likelihood of a star hav- 
ing masses under each side of the measurements' prob- 
ability distribution. This is both because more massive 
stars are intrinsically rarer due to the IMF, and because 
stars of different masses evolve at very different rates. 
The bottom panel of Figure Q] illustrates this difference 
in evolutionary rates for the YREC models with uni- 
form time sampling (At = 1 Myr). Because they evolve 
slower, there are many more grid points for a Solar-mass 
star on the subgiant branch than there are for more mas- 
sive subgiants. 

The prior mass distribution must accurately reflect the 
relative numbers of observable stars of various masses. 
To generate this prior we used the TRILEGAL Galactic 
stellar population simulation code that incorporates the 
IMF, stellar evolution and photom etric system to pro - 
duce synthetic stellar populations ()Girardi et al.l T2005) . 
We accessed the code using Perl scripts provided by L. 
Girardi (2012, private communication), with the default 
Galactic population parameters of the online TRILE- 
GAL 1.5 input form. Given a particular star of interest, 
one can query the simulation to determine the realisti- 
cally expected distribution of masses at a given set of 
observed photometric properties. In our case we use the 
B — V color and absolute magnitude (M.w) since these 
properties are avai lable from the Hippar cos catalog for 
all subgiants in thelJohnson et al.l (|2010al ) target sample 
(jvan Leeuwenll2007| ). 

In Figure [2] we illustrate the effect of incorporating 
such a prior into the mass measurement for a particular 
subgiant. In this case the star has B — V = 1.00 ± 0.02 
and My = 2.2 ±0.5. The stellar mass is M = 1.7 ±0.12 
Mq based on spectroscopy alone, which represents the 
likelihood term. We construct the prior using all stars 
from the TRILEGAL simulations with similar photomet- 
ric properties, using a 1-a cut in Aiw and 3-er cut in 
B — V, to allow for enough stars in the simulation. We 
find that while this prior distribution peaks at 1.15 Mq 



the posterior distribution resulting from the product of 
the likelihood and prior peaks at 1.65 Mq. Including the 
prior, which contains information about the IMF and 
different evolution rates, results in a mass estimate that 
is 3% lower than the likelihood alone. Thus, the differ- 
ent evolution rates of subgiants of various masses is not 
enough to lead to systematically overestimate any indi- 
vidual subgiant's mass by « 50%, as suggested by Lll. 

Note that the prior distribution peaks near 1 Mq , sim- 
ilar to the simulated mass distribution of Lll. This is be- 
cause the prior contains no information about the star's 
actual metallicity. Once a spectroscopic assessment of 
the star's metallicity (and T e g) is made, the likelihood 
function modifies the prior accordingly. That, combined 
with the high precision of the spectroscopic parameters 
results in a likelihood term that dominates over the prior 
and favors higher masses. 

3. THE MASS DISTRIBUTION OF SUBGIANTS IN 
THE KECK DOPPLER SURVEY 

We next turn our attention to the question of whether 
these massive subgiants should exist at all. Lll claimed 
that massive, evolved stars with M > 1.5 M Q should 
be exceedingly rare. So much so that the relative num- 
bers of subgiants with masses in excess of 1.5 Mq com- 
pared to 1.0-1.2 M Q should be taken as evidence that 
the mass measurements must be incorrect. To test this 
hypothesis, we assess the expected distribution of masses 
for stars selected in t he same manner as the targets of 
Uohnson et all (|2010bt ). As we will show, accounting for 
all selection criteria is key to properly estimating the ex- 
pected mass distribution of sample of subgiants. 

Lll simulated stellar populations using the YREC stel- 
lar evolution models together with assumptions about 
the form of the Galactic initial mass function (IMF). 
These features effectively imposed a prior in his Monte 
Carlo simulations as stars were drawn far less frequently 
for masses greater than 1.5 M Q than those closer to So- 
lar. In his simulations, Lll found that only 11% of his 
simulated subgiants had M > 1.5 M Q . 

3.1. The importance of a magnitude limit 

Stellar evolution and the IMF do not have the final say 
in shaping the distribution of stellar masses for Doppler 
survey targets. Surveys of subgiants have specific selec- 
tion criteria that result in stellar samples that are very 
different from the Galaxy's stellar population as a whole. 
For example, the sample of subgiants monitored at Keck 
Observatory were selected based on 0.8 < B — V < 1.1 , 
1.8 < M v < 3.0, and V < 8.5 (| Johnson et al.ll2010al) . 
Another important criterion used to select subgiants is 
the requirement that the s tars have M > 1.3 Mq when 
their Hipparcos B-V colors (jvan L eeuwcn 2007) and abso- 
lute V-band magnitudes (My ) are interpolated onto the 
Padova model grids, under the assumption [Fe/H] = 
(| Johnson et al.ll2010bl ). 

The magnitude criterion was not used in the simula- 
tions of Lll, but it has a profound impact on the ex- 
pected mass distribution of a sample of target stars. This 
is because the flux received from a star scales with dis- 
tance as d~ 2 , while the volume occupied by stars within 
that distance scales as d 3 . Since L ~ M 2 5 on the sub- 
giant branch (based on inspection of the YREC models), 
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subgiants with M — 1.8 M© and visual magnitude V = 8 
observed along a given Galactic sight line will occupy a 
volume [L(M = 1.8)/L(M = 1.2)] 3 / 2 = (1.8/1. 2) 4 « 5 
times larger than a star with M — 1.2 M and the same 
apparent magnitude. This effect offsets the smaller num- 
ber of massive stars in the Galaxy due to the IMF and 
evolution rate. The IMF predicts the number of stars per 
mass interval scales as M~ 2A , and the evolution rate re- 
sults in a relative number of stars on the subgiant branch 
that scales as M -1-8 . Combined, the deficit of massive 
stars in the Galaxy scales as M -4 2 . 

Thus, an apparent magnitude cut will increase the 
number of observable, massive subgiants within a given 
apparent magnitude range, which compensates for the 
dearth of more massive stars due their shorter evolu- 
tion timescale and the stellar IMF. This effect is simi- 
lar to the Malmquist bias in galaxy surveys, in which 
the magnitude-limited survey will result in an apparent 
overabun dance of massive, luminous galaxies at higher 
redshifts (Malmquist 19220- The simple scaling argu- 
ments presented here give a rough sense for the relative 
numbers of stars of various masses within a magnitude- 
limited survey, but they do not account for all effects that 
will ultimately shape the mass distribution. For a more 
thorough analysis we again turn to Galactic population 
models. 

3.2. Simulating the expected mass distribution of 
subgiants 

We estimate the stellar mass prior by first simulating 
samples of stars over the entire sky with a wide range 
of apparent magnitudes. We then select subgiants from 
these simulated samples in the same manner th at the re- 
tired A stars surveyed bv lJohnson et al.1 (|2010b[ ). namely 
0.8 < B - V < 1.1, 1.8 < My < 3, V < 8.5 and the re- 
striction that the stars' colors and absolute magnitudes 
correspond to M > 1.3 M Q based on Solar- metallicity 
stellar models. 

We simulated 768 lines of sight, uniformly distributed 
across the sky using the Hierarchical Equal Area iso- 
Latitude Pixelization (HEALPIX) schem^j, with the ex- 
tinction at infinity c alculated by the NASA/IPA C ex- 
tragalactic database dSchlafiv fc Finkbeinerl 120111 ). We 
avoided the galactic plane (|&| < 4°) because TRILE- 
GAL is known to exhibit "major discrepan cies" with ob- 
serva tional surveys in the Galactic plane (jGirardi et al.1 
120051) , and because of the large number of stars returned 
by those simulations. 

Figure [3] shows the resulting distribution of simulated 
subgiant masses. For ease of future use, we adopt an 
analytic form of the distribution; we find the posterior 
distribution is described well by the expression 

p(M|z)«4=r ex p(-f) W 

V 27ro M V 2 J 

6 As an historical aside, K. G. Malmquist also published a study 
of the dis tribution of (unevolved) A-typ e stars in the Solar neigh- 
borhood (Malmauist & Hufnaecl 1933). In principle, a similar 
study could be used to check the mass measurements of subgiants 
by comparing the ratio of A-type stars to the number of equally 
massive subgiants. This ratio should be equal to the ratio of the 
main-sequence lifetime of A dwarfs to the lifetime of stars on the 
subgiant branch. 

7 http://healpix.jpl.nasa.gov/ 
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Fig. 3. — Simulated sample of stars selected with the same criteria 
as the subgiants in RV surveys (black histogram) . The best-fitting 
log- normal distribution is also shown (red). For comparison, we 
have also shown the mass distribution adapted from Lll (blue 
dashed lines), which is based on a volume-limited survey rather 
than a magnitude-limited survey such as ours. 

where 

log 10 (Af)-Af 

0~ M 

and M = 0.134 M , a M = 0.0783 M and the di- 
mensional skew term is s = —0.253, for all masses 
M > 0.7 M . 

Our stellar mass distribution is qualitatively similar to 
that shown in Figure 4 of Lll. However, while the Lll 
distribution peaks sharply at 1.0 M , ours has a peak 
near 1.3 M , with relatively few Solar- mass subgiants. 
Indeed, there are twice as many stars in our simulations 
with M > 1.5 M as there are for M < 1.1 M . The 
reason for the differences between our simulation and 
Lll's is the V-band magnitude cut, and the a priori selec- 
tion of stars that reside near model tracks corresponding 
to M > 1.3 M Q . While the IMF and subgiant evolu- 
tion rate favor less massive stars, the higher luminosities 
of more massive subgiants makes them visible within a 
much larger volume. 

To test the effects of our added selection criteria com- 
pared to Lll, we selected stars by relaxing certain cuts. 
By relaxing the criterion M > 1.3 M when stars' Aiy 
and B — V colors are compared to Solar- metallicity model 
grids, we find that the low-mass tail of the distribu- 
tion is filled in, which brings the peak of the posterior 
distribution down to 1.2 M . When we impose a vol- 
ume limit, we recover a distribution very similar to the 
volume- limited sample shown in Figure 4 of Lll. 

4. CONCLUSIONS 

We conclude that the findings of I Johnson et alj (|2010aD 
are still broadly correct: there is no compelling reason to 
doubt the existence of a sizeable number of subgiants 
with masses in the range 1.5 M to 2.0 M , and the 
giant planet occurrence rate is much higher around these 
stars than around Sun-like stars. A full recalculation of 
the stellar mass dependence of giant planet occurrence 
is beyond the scope of this work, but will be addressed 
with a larger sample of planet detections out to ~ 4 AU 
in a future contribution. 
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That there are evolved stars with masses significantly 
larger than the Sun's should not be surprising. Within 
only 11 parsecs there resides the slightly-evolved K giant 
f3 Geminorum (=Pollu x), with M ~ 2 Mm, and a long- 
perio d Jovian planet ((Hatzes et al.l 120061 : IReffert et al.l 
2006). While not a subgia nt by the selection crite- 
ria of Uohnson et al.l ()2010bf ). (3 Gem does lie within 
the color and magnitude cu ts defined by Lll, with 
My = 1.08, B — V = 0.99. lHatzes et all ((2013 ) veri- 
fied previous model-based mass estimates by measuring 
the stellar density using photometric and spectroscopic 
asteroseismology. Together with an interferometric an- 
gular diameter and Hipparcos parallax they measured 
M = 1.91 ± 0.09 M . This value compares well to the 
interpolated mass M — 1.96 ± 0.19 M , based on the 
Padova model grids. 

We conclude that within the confines of our current 
best understanding of stellar structure and evolution, 
many of the evolved stars observed in Doppler surveys 
are indeed as massive as A stars. There remains the pos- 
sibility that there exist systematic errors in the stellar- 
evolution models that are not accounted for in our tests. 
However, these systematic errors cannot be too severe. 
Otherwise, it would be difficult to explain the close match 
between the empirical and model-based mass measure- 
ments of /3 Gem. It would also be difficult to explain 



how lSandage et all (|2003[ ) were able to measure ages of 
the Galaxy's oldest stars by fitting isochrones to the lo- 
cal subgiant branch and find an answer consistent with 
the age of dwarf stars near the local main-sequence turn- 
off. Also, as noted by Lll, it would be especially difficult 
to explain the apparent correlation between stellar mass 
and the occurrence of giant planets. 

Nevertheless, tests of systematic errors in stellar evo- 
lution models using planet transit light curves, eclips- 
ing binaries and asteroseismology are very much worth- 
while. Fortunately, the large number of transiting plan- 
ets and ecli psing binaries in the NASA Kepler mission 
target field (jPrsa et aT1l2011l ). together with the exquisite 
photometric precision produced by the Kepler space tele- 
scope, will provide many opportunities for these tests in 
the near future. 
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